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ABSTRACT 

The presence of a close, low-mass companion is thought to play a substantial and per- 
haps necessary role in shaping post- Asymptotic Giant Branch and Planetary Nebula 
outflows. During post-main-sequence evolution, radial expansion of the primary star, 
accompanied by intense winds, can significantly alter the binary orbit via tidal dis- 
sipation and mass loss. To investigate this, we couple stellar evolution models (from 
the zero-age main-sequence through the end of the post-main sequence) to a tidal 
evolution code. The binary's fate is determined by the initial masses of the primary 
and the companion, the initial orbit (taken to be circular), and the Reimers mass- loss 
parameter. For a range of these parameters, we determine whether the orbit expands 
due to mass loss or decays due to tidal torques. Where a common envelope (CE) phase 
ensues, we estimate the final orbital separation based on the energy required to unbind 
the envelope. These calculations predict period gaps for planetary and brown dwarf 
companions to white dwarfs. The upper end of the gap is the shortest period at which a 
CE phase is avoided. The lower end is the longest period at which companions survive 
their CE phase. For binary systems with 1 Mq progenitors, we predict no Jupiter- 
mass companions with periods ^270 days. Once engulfed, Jupiter-mass companions 
do not survive a CE phase. For binary systems consisting of a 1 Mq progenitor with 
a companion 10 times the mass of Jupiter, we predict a period gap between ~0.1 and 
~380 days. These results are consistent with both the detection of a ~50 Mj brown 
dwarf in a -0.003 AU (--0.08 day) orbit around the white dwarf WD 0137-349 and 
the tentative detection of a —2 Mj planet in a >2.7 AU (>4 year) orbit around the 
white dwarf GD66. 

Key words: stars: AGB and post-AGB - stars: low-mass, brown dwarfs, planetary 
nebulae: general 



1 INTRODUCTION 

For low- and intermediate-mass stars (initially <8 Mq), 
post-main sequence (post-MS) evolution is characterized 
by expansion via giant phases accompanied by the onset 
of mass-loss. In particular, during the Asymptotic Giant 
Branch phase (AGB), dust-driven winds expel the stellar 
envelope as the star begins its transition to a white dwarf 
(WD). However, before formation of the WD remnant, the 
spherical outflows observed during the AGB phase undergo 
a dramatic transition to the highly asymmetric and often 
bipolar geometries seen in the post-AGB and planetary neb- 
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ula phases (PN; Sahai and Trauger|1998 1. This transition is 
often accompanied by high-speed, coUimated outflows. For 
recent reviews se e [Balick and Frank] ( |2002[ ), [van Winckel| 
( 20031), Ide Marco|(|2009|). 



A central hypothesis to explain shaping in post- 
AGB/PNe is that a close companion is necessary to power 
and shape bipolarity. This is supported by observations of 
excess momenta in almost all post-AGB outflows relative 



to what isotropic radiation pressure can provide (Bujarra- 
|bal et aL]|2001[ ). Additionally, close companions have been 
detected in or around giants, bipolar post-AGB and PNe 
systems and in some cases seem to be responsible for out- 



flow shaping ( 


Silvotti et al.| 


2007 


De Marco et al.|2008 


Sato| 


let al.||2008b|a| |Witt et al. 


12009 


IMiszalski et al.||2009b|a| 
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Niedzielski et al. 2009 Chesneau et al. 2009 1 . In particu- 



lar, the recent detection of a white dwarf with an orbiting 
~50 Mj brown dwarf in a ~2 hour orbit demonstrates that 
low-mass companions can survive a common envelope phase 
(CEP) ( [Maxted et~aI1|2006[ ). The detection of a planetary 
companion (Msini = 3.2 Mj) around the extreme horizon- 



tal branch star V391 Pegasi in a 
period 



-1.7 AU orbit (~3.2 year 



Silvotti et al. 20071 and the tentative detection of 



a ~2 Mj planet in a >2.7 AU orbit (>4 year period; MuL 



lally et al.|2008||2009| ) around the white dwarf GD66 provide 
motivation for this work. 

Indirect observational evidence for binarity comes from 
maser observations that suggest magnetic jet coUimation in 



AGB and young post-AGB stars (Vlemmings et al. 2006 



|Sabin et al.||2007| [Vlemmings and van Langevelde||2008| ) 

Such collimation supports the binary hypothesis because it 
is difficult for single AGB stars to generate the large field 
strengths necessary to power the outflows without an addi- 
tional source of angular momentum ( Nordhaus et al.]|2007 
2008a I. If a close companion is present, strong interactions 



can transfer energy and angular momentum from the com- 
panion to the primary. In particular, if the companion is 
engulfed in a common envelope (CE), rapid in-spiral can 



cause significant differential rotation of the envelope ( Nord- 
haus and Blackman||2006| [Nordhaus et al.||2008b[ ). Coupled 



with a strong convective envelope, large-scale magnetic fields 
are amplified and are sufficient to unbind the envelope and 
power the outflow (Nordhaus et al.|2007|. 



Moe and De Marco 



( 2006 1 overpredicted the Galactic 
PN population by a factor of ~6 for PN with radii <0.9 pc 
(discrepant at the 3 a level; Jacoby[[l980 Peimbert][1990 1. 
This discrepancy could be alleviated if only a fraction of 
stars become PNe. The authors argue that perhaps inter- 
acting binaries or, in particular, common envelope systems 
may be responsible for producing the majority of Galactic 
PNe ( De Marco and Moe|2005 L However, to determine the 
validity of a hypothetical CE/PN connection, it is impor- 
tant to understand which binary systems will undergo a CE 
phase during their evolution. The fact that low-mass stel- 
lar companions can survive CE evolution motivates further 
study independent of a potential connection to PNe ( [Maxted] 
et al.|2006| |. 

Various aspects of tides on the evolution of post-MS bi- 



naries have previously been investigated (Rasio et al.[[1996 



Carlberg 


et al. 


[2009] 


Carlberg 


et al.[ 


12009^ 



Villaver and Livio 2009 1 . Recently, 



aimed to understand fast rotation in 



trying to predict the distribution of planets around evolved 
stars. [Carlberg et al.[ ( [2009[ ) calculated orbital evolution sce- 
narios but neglected mass-loss effects. [Villaver and Livio[ 



(20091 included stellar mass-loss in their model and mod- 
eled a range of progenitor masses (from 1 to 5 Mq), but 
did not examine how results depend on the prescriptions for 
either tidal interactions or stellar mass loss. They consid- 
ered a wide range of possible astrophysical effects, including 
frictional and gravitational drag, wind accretion, and atmo- 
spheric evaporation. We neglect these effects in our study; 
most are negligible compared to tidal torques and mass loss 
from the primary. The exception is perhaps evaporation, 
which has been claimed to destroy planets of mass < 15Mj 
( [Villaver and LiviopOOT] [Livio and Soker|1984| ). This is sim- 
ilar to the minimum mass required to unbind the stellar 



Q. 
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Figure 1. The period gap for low-mass companions around white 
dwarfs. The orbit of a companion located initially at ai decays 
and plunges into the giant star. Depending on the mass of the 
companion and stellar structure at plunge time, the companion 
may or may not survive the CE phase. Companions slightly ex- 
terior to ai avoid engulfment and never enter a CE; their orbits 
expand due to mass-loss. The gap is set by the final maximum 
semimajor axis which survives CE evolution (a,nin) and the final 
minimum semimajor axis which avoids tidal engulfment (amax)- 



envelope (Section 6). Our work builds upon both previous 
results in that we follow the evolution from the zero-age 
main sequence (ZAMS) through the entire post-MS, and we 
consider a variety of prescriptions for tidal torque and for 
stellar mass-loss rates. 

In this paper, we investigate how low-mass companions 
(<0.1 Mq) become immersed in a common envelope. By 
coupling stellar evolution models to tidal evolution equa- 
tions, we can determine when (i) mass-loss dominates and 
leads to orbital expansion or (ii) dissipation of tidal energy 
dominates and leads to orbital reduction. In particular, for 
a given orbital separation, we can determine which binary 
systems incur a CE by plunging into their host stars modulo 
the uncertainties in the theories of tidal dissipation and stel- 
lar mass-loss. In ^ we describe our equations and discuss 
a possible tidal dependence on period. In !|3j we present our 
stellar evolution models and focus on low-mass companions 
(planets, brown dwarfs and low- mass main sequence stars). 
In ^ we present our results and discuss the implications 
of two commonly used tidal prescriptions. For systems that 
incur a CE, we use an energetics argument to determine the 
maximum radius at which a given companion can survive the 
interaction by successfully ejecting the envelope. Further- 
more, we determine the minimum separation that evades a 
CE phase. By explicitly following the orbital evolution, we 
calculate the final separation of the binary system. For a 
given binary, this produces a gap inside of which we expect 
the absence of companions. From the set of all gaps, we can 
determine a minimum period gap for planetary companions 
around white dwarfs (see Fig. [T]). The prediction of a gap 
should be of interest to recent searches for both substellar 
and, in particular, planetary companions to white dwarfs in 
addition to future observations ([tluUally et a"L][2008[ [2009] 



Kilic et al.|2010 Qian et al.|2010 1. We discuss these results 
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in ^ and !|6] In ij7j we discuss the possibility of detecting 
WD+companion transits and applications to GALEX data. 
We conclude in ^ 



2 TWO-BODY TIDES 

We employ a two-body gravitational and tidal interaction 
model that consistently couples the evolution of the mass 
and the radius of the primary star with the orbit of the 
companion. It includes the tides raised on the companion (by 
the star) and the tides raised on the star (by the companion). 

The evolution of the semimajor axis (a) is due to two 
contributions: the tidal interaction between the primary and 
the companion and the mass loss of the primary 



da 



^ tides 



(1) 



We employ the equations of Ferraz-Mello et al. ( 2008 1, which 



use the Q-formalism of Goldreich ( 1963 1 for the tidal evo- 



lution of the orbital eccentricity, e and the semimajor axis, 
a. This approach has been used extensively to model the in- 
flated radii and orbits of transiting extra-solar giant planets 



(EGP; 


Kaula||1968 


Ferraz-Mello et al.||2008 


Jackson et al. 


2008 


I 


bgui and Burrows|2009||Miller et al.|2009||Ibgui et al. 


2009 


Jackson et al.|2009l. 



Our formalism assumes that the primary star and com- 
panion are both spherical and that their rotational and or- 
bital angular momenta are aligned. Furthermore, we assume 
the companion spin is synchronized with its orbital period 
and that the mean orbital motion is greater than the stel- 
lar rotation rate. This is reasonable as typical equatorial 
surface velocities are ~1 — 3 km s~^ for low-mass giants im- 
plying that, for companions that are close enough for tides 
to matter in the post-main-sequence evolution, the angular 
velocity of the orbit is much greater than the angular veloc- 
ity of the primary star (Gray||1989j [Massarotti et al.|2008 |. 



Under these assumptions, we have that 



1 da 


1 


a dt 


ai3/2 


1 de 


1 


e dt 


ai3/2 



_8_ 
25 



57 2 



(2) 



where Rc and 7?* are the companion and primary radii, Q'^ 
and Q'^ are the tidal dissipation factors of the companion 
and the primary and A'1,2 are constants given by 



K2 



63 



\l/2 



Me 



\l/2 



4 

225 

where Mc is the companion mass and K'h 



(3) 

is the primary 



Furthermore, the radius and mass of the companion are 
held constant while those of the primary are time-dependent. 
Note that in these equations, e can only decrease and if the 
primary mass were constant, a could only decrease as well. 

The tidal dissipation factors Qj, and Q'^ are dimension- 
less parameters given by Q' = 3Q/2fc2, where Q is the spe- 
cific tidal dissipation function and k2 is the Love number 
(|Goldreich|1963| [Goldreich and Soter|1966[[Ogilvie and Lin| 
20071. Q is a dimensionless parameter that quantifies the 



efficiency of tidal dissipation inside a body. It is defined by 
its reciprocal, the specific dissipation function 



2-KEa 



dE 
~dt 



dt, 



(4) 



where Eo is the maximum energy stored in the tidal dis- 
tortion of the body and the integral ,f dt is the 
energy dissipated over one orbital cycle (Goldreich 19631. 
Typical estimates yield Q'^ ~ 10'^ °"*-°, Q'^ ~ ;lqo.o-d 5 
for Jupiter mass pla nets around main-sequence stars and 



( Goldreich and Soter 



1966 



Qjupitc . 

,Poalc^l98l||Ogilvie and Lin|2004||2007||Jackson"et al..2009| 



Yoder and 



2.1 Primary and Companion Tides 

Tides on the primary act to synchronize its spin while tides 
on the companion act to circularize the orbit. We can com- 
pare the tides raised on the companion (induced by the pri- 
mary) to the tides raised on the primary (induced by the 
companion) as follows: 



jda/dt)^ 
{da/dt)^ 

jde/dt)^ 
(de/dt)^ 



7e^ 



28 
25 



fe2 



Q'c 



(5) 



(6) 



For a 1 Mq AGB star (with radius ~1 AU) and a 0.1 Mq 
main sequence companion (with radius ~0.1 -Rq), with e — 
0.1, we have 



(7) 
(8) 



For lower mass companions, these ratios remain extreme. 
Therefore, for Q'^ ~ Q'^, if the primary is a post-main se- 
quence star, the orbit decays almost entirely due to tides 
raised on the star, not to those raised on the companion. 

As the orbit circularizes (e — >■ 0), the tides on the com- 
panion decrease until they vanish for perfectly circular or- 
bits. However, note that some dissipation must occur within 
the companion to maintain synchronous rotation. In this 
limit, the orbital evolution is governed by non-synchronous 
rotation between the companions orbit and spin of the pri- 
mary. As a result, in-spiral continues and might lead to the 
companion plunging into the primary (the cases we are in- 
terested in here). 

2.2 Main-sequence and Post-MS Tides 

For the known extrasolar giant planets (EGPs), the tides 
raised on the primary can be substantial and can lead to sig- 



nificant orbital reduction during the main-sequence (Ibgui 
and Burro^^|200^ 



Levrard et aL] [20091 [Ibgui et al] 



2009 



Spiegel et al.[ 20101. For the same Qi, we can estimate 
whether a similar reduction is expected during the post- 
MS. From the previous section, we have that the tides on 
the primary dominate and, hence. 



Alna 



rgb 



AlnOn 



Rr 



Rms 

Mc,rgb 
Mc.ms 



-13/2 



(9) 



*,rgb 
k.agb 
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Figure 2. Reproduced versions of Fig. 4c of | |Verbunt and Phin-| 
|ney||1995l l assuming a period dependence on /. The transition 
between circularized and non-circularized systems is preserved 
around unity for a variety of scalings. Here we present one exam- 
ple. Left: / oc which implies Q'^ oc P. Right: / oc which 
implies Q'^ oc P". 



where Tms is the typical lifetime during the main sequence 
and Trgb is the lifetime during the RGB phase. 

For the main sequence, we take a 1 Jupiter mass com- 
panion around a 1 Mq main-sequence primary with semi- 
major axis ~0.05 AU. For the RGB phase, we assume that 
the companion is orbiting at 1.3 AU, that Q^t.^s ~ Q*,rgb 
and that -Rrgb ~ W^Rms and Trgb ~ lO^^rms. Using these 
parameters yields AlnOrgb/ Alnams ~ 0.063. Therefore, for 
this system, we expect tides to be weaker on the RGB phase 
than they are for EGP systems. This is likely to be true for 
low-mass primaries (< 1.5Mq) where radial expansion dur- 
ing the RGB is substantial (see Fig. [sj. For higher mass 
primaries, radial expansion is minimal during the RGB and 
extensive during the AGB. However, typical AGB lifetimes 
are ~10~^ — 10~''Trgb. We can estimate the effect of tides dur- 
ing the AGB phase in the same way. Assuming i?agb ~ 200 
^0, Q*,ms = Ql.agb and Tagb ~ lO'^Tms, we havc that 
AlnOagb/ AlnOnis ~ 2.0 X 10"'^. Thus in general, we expect 
synchronization tides on the AGB to be weaker both than 
those on the RGB and for typical EGPs around main se- 
quence stars. Note the steep dependence on the ratio of the 
stellar radii in Eq.jo] If -Ragb = 2i?rgb and Tagb ~ (l/30)rrgb, 
then Alnttagb/ AlnOrgb ~ 1 and tides are approximately 
equal during the RGB and AGB phases. However, as will 
be seen shortly, there is reason to expect that may be 
orders of magnitude smaller on the giant branches than on 
the main sequence for the same a/i?*, with a correspond- 
ing increase in the importance of the advanced evolutionary 
phases for the orbit. 



2.3 A Period Dependence of i 



It is worth nothing that the work of Goldreich ( 1963 1 and 
[Goldreich and Soter| (|1966[) doe s not specify a tidal dissi- 
pation mechanism. Zahn ( 1966 1 however, proposed a tidal 



theory based on turbulent viscosity. This theory was tested 
and calibrated in stellar binaries that contain an evolved star 
primary (Verbunt and Phinney|1995 1. Evolved stars are ex- 
pected to possess extended convective zones, which could be 
crucial for tidal dissipation. In their formalism, the authors 
introduce a dimensionless factor / that is calibrated via the 
orbital eccentricity measurements in their post-MS binary 
sample. The authors argue that observational data imply 
that / is constant and approximately unity. 

By setting equal the expression for the tidal torque (as 
a function of Qi) in Goldreich and Soter ( 1966 1 to the ex- 
pression (as a function of /) in Zahn (19891, as simplified 
by Verbunt and Phinney ( 1995 I, we have that 



01 



63 GM* 



16tt 

'Tconv 
'dyn 



P 



A'/cnv -R* 

L 



1/3 



(10) 



where Mcnv is the mass of the convective envelope, P is the 
orbital period, Tconv is the convective timescale and rdyn is 
the dynamical timescale. If / is in fact constant, then Q'^ 
is proportional to orbital period. Note that the numerical 
prefactor in Eq. (10) (63/167r) is slightly different from the 
225/327r that results from setting our expresion for dhie/dt 
(in Eq. 2) equal to the corresponding expresion in Verbunt 



and Phinney (19951. The difference (compared to the vari- 



ation in period) is sufficiently small that it does not effect 
our results. 

The argument for / constant and ~1 derives from 
Fig. 4c of [Verbunt and Phinney] ( |1995[ ). The transi- 
tion between circularized and non-circularized systems 
occurs when the abscissa is 0. The abscissa is given 
by logio{-Aln[e]/f} = logio {In [ci/et]} - logio {f}. If 
logic {In [ci/ef]} is greater than (less than) 0, the sys- 
tem is strongly (barely) circularized. The transition oc- 
curs near logio {In ei/ef]} = 0. Looking at Fig. 4c of |Ver-| 
bunt and Phinney] ( 1995 1, we see this transition happens at 
logio {In [ci/ef]} - logio {f} ~ 0, therefore, / ~ 1. 

While this argument does constrain / for systems 
at the sharp transition between circularization and non- 
circularization, it does not rule out a dependence of / on 
period. In fact, the transition between circularized and non- 
circularized systems remains at log {A In [e]} ~ if / scales 
as P^~^ with 1 > a; > and has an appropriate normaliza- 
tion (see Fig. [2|. This corresponds to Q'^ oc P^ . It is worth 
mentioning that the scaling used in Fig. [2] is one example 
among many that preserve the location of the transition. 
Different normalizations preserve the transition for different 
ranges of x. For example, / = (P/200days)^~^, preserves 
the transition for 2 > x > 0. T he value x = 2 (for which 
Ol oc P~^) has been urged by Goldreich and Nicholson 
( [1977 1 in the limit P <^ TTTconv; however, P > Trrconv for 
most of the cases of interest to us. 

Thus, by using Eq. (|10[), we are able to employ the 



(pei 



tidal prescription of Zahn (11966 1 with the observational cal- 
ibration of Verbunt and Phinney (19951. In Section 4, we 



show that employing such a formalism leads to values of Ql 
between ~10^ and 10^ for post-main sequence giants. These 
values are two to seven orders of magnitude lower than those 
typically invoked in the context of extra-solar giant planets 



(Goldreich and Soter 1966 Yoder and Peale 1981 Ogilvie 
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and Lin||2004l [2007| | Jackson et alT] [20091 [Tbgui and Bur-| 
rows|2009| [Ibgui et al.||2010| [Miller et al.||2009[ ) and lead to 
strong tidal interactions. We pause to note that tidal theory 
is an active field of research (Goodman and Lackncr 2009[ 



Gu and Ogilvie|2009 Arras and Socrates|2009 1 Particularly 
relevant is a recent study that investigated a frequency de- 
pendence of Q' ( Greenberg[[2009 1. This work demonstrated 
that commonly used analytic approaches (which make anal- 
ogy to a driven harmonic oscillator) are probably valid only 
for low eccentricities and inclinations. Theoretical and ob- 
servational constraints on the tidal dissipation mechanism 
will help constrain the results of this work in the future. 

For now, in light of the uncertainties in Q'^, and in /, 
we adopt the following form 



X 



(11) 



where 11 = P/r^y-n = P (pG) is the orbital period divided 
by the dynamical time, Ho is a reference value, logjQ Qo is 
between 5 and 9. For Ho, we use a 10 day orbital period 
divided by the dynamical time of the Sun, namely Ho — 
(10 days) x (pqG)^^^ . Note that even in the case where f — 1 
and X = 1, Eq. (Ill is not fully consistent with Eq. ( 10 1. For 



this paper, we focus on two representative values of x such 
that oc n° and oc n\ 



3 PRIMARY AND COMPANION MODELS 

Our stellar evolution models are calculated using the 
"Evolve Zero-age Main Sequence (EZ) code" ( Paxton|2004 |. 
For each initial stellar mass, we evolve the star from the 
ZAMS through the end of the AGB phase including a 
Reimers mass-loss prescription fReimers"1975 1. In order to 



explore the influence of the rate of mass-loss on our calcu- 
lations, we use a range of values of the Reimers rj param- 
eter: 0.7, 1, and 5. Each stellar model has a metallicity of 
Z = 0.02. In addition to the stellar radius and mass, we cal- 
culate the core mass and envelope binding energy as a func- 
tion of time. Our stellar models (for rj — 1) are presented in 
Fig. [3| while key parameters for all models are summarized 
in Table [H 



3.1 Companions 

We consider three types of companions: planets (Mc < 



O.OO26M0; Zapolsky and Salpeter 19691, brown dwarfs 
(O.OO26M0 < Mc < O.O77M0; [Burrows et al.[[1993[ ) and 
low-mass main sequence stars (0.077 < Mc < O.IMq). We 
adopt the following approximation to the models of |Burrows| 
et al.[ ([1993| [1997[ i2001) for brown dwarf radii 



R2 = ( 0.117- 0.0541ogio 

Mc 

0.0026 



Mc 



-|-0.0241ogi 



0.0026 

Rq , 



(12) 



and use a homologous power-law for low-mass stellar radii 
R.-l^) Re, (13) 
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Figure 3. Mass (dotted line) and radius (solid line) profiles for 
a subset of our stellar models (i.e. those with t] = 1). The mass 
profiles have been offset 5% in time so that the evolution is dis- 
tinguishable on the plot. 



We carry out our calculations for the following compan- 
ions: a 0.001 Mq (1 M,j) planet, a 0.01 Mq (10 Mj) object 
that could be either a massive planet or a low-mass brown- 
dwarf (depending on formation scenario; hereafter referred 
to as a brown dwarf), and a 0.1 Mq (100 Mj) low-mass main 
sequence star. 



4 TIDAL RESULTS 

As the binary system evolves, mass-loss and tidal torques are 
in competition. Mass lost from the system acts to increase 
the semimajor axis while tidal torques decrease it. For each 
primary and companion, we compute the evolution of the 
orbit from the zero-age main sequence through the post- 
main sequence. If the companion is tidally captured (i.e. 
plunges into the primary star) , it enters a common envelope 
with the primary and we halt the calculation of the orbital 
evolution. If the companion evades tidal engulfment, mass- 
loss continues and the orbit expands until the end of the 
evolutionary model. 

We performed calculations for two tidal prescriptions: 
those corresponding to Eq. ([2| and those corresponding to 
Eq. (10 1 when / = 1. With this approach, we can directly 
compare the tidal theories of [Goldreich and Soter[ ( [1966[ ) to 
those of Zahn ( 1966 1 when / is constant and equal to one 
(Vcrbunt and Phinney|1995[ ). 

For each stellar model, companion mass, and tidal 
theory, we calculate the maximum initial semimajor axis, 
aii,max, that is tidally engulfed. Eq. ^ is accurate to low- 
est order in eccentricity and valid for e <C 1 ( [Goldreich and[ 
|Sotcr 1966). For low e (<0.1), the effect of eccentricity is 
to modify the dominant tidal term by <10%. Thus, we as- 
sume circular orbits for the ensuing calculations. It should 
be noted that the formalism we employ only considers the 
quadrupole component of the tidal potential. 

Although the Reimers prescription is an imperfect de- 
scription of mass-loss, we investigate how different rates of 
mass-loss influence our results by using rj values of 0.7, 1, 
and 5 ( [Reimers[1975[ ). Our mass-loss prescription reproduces 
typical mass-loss rates during the RGB but underestimates 
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3.00 


0.20 


2.99 


3.71 


1.60 


0.43 


1.07 






1.00 


0.87 


0.73 


0.75 


0.56 


0.46 


0.56 






1.50 


0.67 


1.38 


1.85 


0.69 


0.45 


0.69 


V = 


1 


2.00 


0.27 


1.97 


2.57 


0.80 


0.39 


0.81 






2.50 


0.16 


2.49 


3.14 


1.07 


0.38 


0.91 






3.00 


0.20 


2.98 


3.50 


1.40 


0.44 


0.99 






1.0(P] 


0.41 


0.38 






0.37 




V = 


5 


2.0() 


0.30 


1.84 


1.37 


0.66 


0.39 


0.60 






3.00 


0.20 


2.94 


2.50 


0.81 


0.43 


0.75 



" Maximum radius on the RGB/AGB. 

^ Calculated when the stellar radius is largest in the RGB and AGB phases. 

The core masses are defined by composition rather than degeneracy. Explicitly, it is the mass interior to 
the outermost hydrogen-free shell. 

The ?7 = 5, 1 Mq model does not achieve a sufficient core mass to undergo a helium flash. 

Table 1. Stellar model data. Column quantities from left to right: zero-age main sequence mass, maximum radius on the RGB, stellar 
mass at the tip of the RGB, maximum radius on the AGB, stellar mass at the tip of the AGB, core mass for the corresponding -Rrgb, 
core mass for the corresponding Hagb. 



mass-loss rates on the AGB. As such, the results presented 
here serve as an approximate upper limit to the parameter 
space where tidal torques dominate over mass loss. Access 
to accurate stellar evolution models and realistic mass-loss 
prescriptions (motivated by observations during all phases of 



stellar evolution) will refine this work in the future ( Schroder 
and Cuntz|2005l[2007t . 



4.1 When / = 1 



The tidal theory of Zahn ( 1966 1 under the assumption that 
/ = 1 has been employed to determine the effect of tidal 



torques on various aspects of post-MS evolution ( 


Soker 1995| 


1996; Villaver and Livio|2009||Carlberg et al.|2009; Boar and| 


Soker 2009 


I. This is accompanied by results that suggest 



strong tidal interactions for companions within ~(5 — 7) x 
^?max, where 7?max is the maximum radial extent of the star 
during either the RGB or AGB phases (|Soker||1995| |1996[ 
pebes and Sigurdss on 2002' lMoe and De Marco '2006 ) . Here, 
we make the same assumption and determine the maximum 
orbital separation that plunges into the primary star for 
various primary and companion masses. 

Figure|4]shows values at the time of tidal engulfment 
under the assumption that / = 1. Typically, values range 
from 10^ "^ < Q'i, < 10^'^, and, for a given progenitor mass, 
Q'i, values vary by less than a factor of ~20, irrespective 
of the choice of ri and companion mass. This is striking, 
as it implies that the assumption of / = 1 leads to tidal 
dissipation rates that are ~2-7 orders of magnitude larger, 
in a dimensionless sense, than those acting on extra-solar 



are typically taken to range from 10^ — 10* ( 


Barnes et al. 


20081 IJackson et al.||2009; Ibgui and Burrows 2009| |Ibgui 


et al.|2010 


Miller et al.|2009 


Levrard et al.|2009 


I 



With small Qf, values, one would expect tidal engulf- 
ment to occur at large separations. Figure [5] shows the max- 
imum semimajor axis that is tidally captured (ai.max) for a 



1000?^^ 



Q'. at plunge-time (f=1) 

M 

c 

M 

c 

M 




Figure 4. Q'^ values for post-main sequence binaries if / = 1. 
Note that choice of / = 1 yields Q'^ values 2 to 7 orders of mag- 
nitude lower than typical values used to explain the orbits and 
radii of transiting extra solar planets. 



range of primary models and companion masses. The top 
panel presents Oi^max in units of AU for the case of / = 1, 
while the bottom panel shows the same results for the case 
of — 10^ (see Table 1). Generally, the more massive stel- 
lar models extend to larger radii on the AGB and, therefore, 
can swallow companions that are farther away. Jupiter-mass 
companions within <2 x iimax are captured while more mas- 
sive companions can plunge at farther distances (10- Jupiter- 
mass companions are captured within a < 2.5 x 7?max; 100- 
Jupiter-mass companions are captured within a < 3 x -Rmax) . 
It should be noted that the time of tidal engulfment roughly 
corresponds to the time at which the primary star is at its 
maximal radial extent. However, this is not exact as most 
binary systems are engulfed when the primary is still ex- 
panding, but near maximum radial extent. 

If tides acting during the post-main sequence are char- 
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■' Maximum a. that plunges (Q .=1 0^ 




Figure 5. Maximum semimajor axis, Oi max that is tidally en- 
gulfed under the assumption that / = 1 (top) and Q'^ = 10^ (bot- 
tom). Both figures present results for companions of mass IMj 
(blue curves), lOAfj (purple curves) and lOOJV/j (red curves). Re- 
sults for different mass-loss models are shown by r; = 0.7 (dashed 
curves), rj = I (solid curves), r/ = 5 (dot-dashed curves). 



duces larger ai^max than the constant case. Finally, for 
the EGPs, it may be that a; < 0, consistent with the pre- 
diction that X should switch from +1 to —1 at P <^ Tconv 
( Goldreich and Nicholson|1977 |. 

In order to estimate the effects of mass- loss, Fig.[7|shows 
ffli.max for T] = 0.7 (left figure) and ?7 = 5 (right figure). In 
general, less mass loss leads to tidal engulfment at slightly 
larger distances. 

By coupling stellar evolution models with various tidal 
theories, we determined the maximum separation at which a 
companion might be tidally captured. Companions slightly 
exterior to ai,max are never engulfed by the envelope of the 
primary but their orbital dynamics are still subject to the 
effects of mass loss and tidal torques. By continuing to follow 
the orbital evolution, we can predict the minimum final sep- 
aration for each binary configuration. The minimum of this 
set is the minimum separation exterior to which one would 
expect to find planetary companions around white dwarfs. 

In the following two sections, we determine the maxi- 
mum separation at which a companion can survive CE evo- 
lution. Interior to this separation, we would expect to find 
companions to white dwarfs. Taken together, the minimum 
and maximum form a separation gap between which we ex- 
pect an absence of planetary and brown dwarf companions 
around white dwarfs. 



4.3 What About the Solar System? 

For the solar system, tides and mass-loss compete and de- 
termine whether the inner planets plunge into the Sun or 
evade tidal engulfment ( Rybicki and Denis|2001 Sackmann 
|et al.|[T993[ ). Under the assumption / = 1, Venus plunges 
into the Sun for all of our 1 Mq evolutionary models. Earth 
however, evades tidal engulfment with the orbit expanding 
for all of our evolutionary models, consistent with the re- 
sults of Rasio et al. ( 1996 1. In addition, none of the planets 



beyond Earth are swallowed. 



acterized by similar Q' values (in the Goldreich and Soter 
1966 formalism) to those acting on EGPs during the main 
sequence (i.e., > 10^, Q'c ~ 10®), then results might be 
significantly different from those derived from the formalism 
of [Zahn| ( [19661 ) ^i^li / = 1 (in Fig- [sf • To compare, in the 
next section we present cases in which Q'^ is constant and 
Q' oc n. 



4.2 Ql constant and oc 11 

The left columns of Fig. [6] show ai^max as a function of 
primary mass for constant Q'^. The right columns present 
cii.max for Q'^, = Qo X (fl/IIo). The evolutionary models in 
Fig. [6] were calculated with a Reimers mass- loss parameter 
of r; = 1. The left plot is in units of AU while the right 
plot is in units of i?max. As such, we can directly compare 
the results of a constant Q'^ to those of / = 1 shown in 
Fig. [5] Note that for 10^ < Q* < 10^ tidal torques are much 
weaker and result in smaller ai,max. In addition, the ratio 
ai,max/^imax < 1 for most of the parameter space. For the 
same companion and stellar evolution model, oc 11 pro- 



5 COMMON ENVELOPE EVOLUTION 

Upon tidal engulfment, the companion enters a common 
envelope w ith the primary star ([Paczynski||1976| |Iben and] 
Livio|1993[|Nordhaus and Blacknian|2006[ ). The velocity dif- 
ference between the orbital motion of the companion and 
the common envelope generates drag. The resulting loss of 
orbital energy leads to rapid in-spiral on week- to month- 
timescales ( Nordhaus and Blackmanl|2006 1. The orbital en- 
ergy released as the companion in-spirals can be used to 
overcome the binding energy of the envelope. If sufficient 
energy is released from the orbit during in-spiral, the com- 
panion can eject the envelope and survive the CE. This is 
expressed as the following: 



-Ebind = —aAEo 



(14) 



where iJbind is the binding energy of the envelope, a is the 
fraction of orbital energy that goes toward ejecting the en- 
velope and AEorb is the change in orbital energy of the 
companion. 

The change in orbital energy of the companion is given 
as -Eorb.j?* — Sorb, a, where Borb.fl* is the orbital energy at 
tidal engulfment (a = i?*) and Eoib,a is the orbital energy at 
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semimajor axis a inside the CE (i.e. a < Rt). The orbital en- 
ergy at tidal engulfment is given by the sum of gravitational 
and potential energies, namely: 



2R-, 



(15) 



where M*[_R*] is the total mass of the primary. The orbital 
energy at a < J?* is given by the kinetic energy [Ka) plus 
the gravitational potential energy {Ua)- 



-Eorb,a ~ Ka + Ua 

GM^ [a] Mc GM^ [R^] Mc 



(16) 



2a 



R* 



^* GM* [r] Mc 



dr 



where the integral term is the work against gravity required 
to move the companion from a toii*. The total change in 
orbital energy is then given as: 



(17) 



GM4a]Mc _ GM4R^,]Mc _ GM-,[r]Mc 
~ 2~a 2R. A 

The companion in-spirals until it is either tidally 
disrupted or supplies enough orbital energy such that 
-aA^Boib = £^bind ( [Nordhaus and Blackmaril [2006] ) . The 
binding energy of the stellar envelope is computed by sum- 
ming the gravitational and thermal energies, at each evo- 
lutionary timestep. We require that the companion supplies 
enough orbital energy to overcome the binding energy of the 
entire envelope. If the companion only ejected the mass exte- 
rior to its orbit, the interior mass would expand to re-engulf 



the companion and restart in-spiral. Note that Eq. (171 de- 



pends on the stellar structure at the onset of tidal engulf- 
ment. 



6 PERIOD GAPS FOR PLANETS AND 
BROWN DWARFS AROUND WHITE 
DWARFS 

To calculate a minimum period gap expected for a given bi- 
nary system, we assume a = 1. This gives an upper bound 
on the orbital radius at which we would expect to find com- 
panions which have survived a common envelope phase (see 
Fig.[l|. In §4, we determined the time at which a companion 
of mass Mc is engulfed by the giant star. Coupled with the 
structure of the star at the time of tidal engulfment, we can 
determine which semimajor axis a satisfies Eq. (14 1. If the 



companion avoids tidal disruption, then it has successfully 
ejected the envelope and survived the CE phase. 

The tidal shredding radius can be estimated by bal- 
ancing the differential gravitational force across the com- 
panion with its self gravity. This yields a tidal shredding 
radius given by Os — Rc \/2M /Mc where Rc is the radius of 
the companion. For a IMj companion around a proto-white 
dwarf core, ~ 7.4 x 10^" cm. For a 10 Mj companioij^ 
Os > 3 X 10^" cm with the actual values dependent on the 
degenerate core mass during the CEP. If the companion un- 
binds the envelope exterior to as, then we say it has survived 
common envelope evolution. Note that the factor of 2M/Mc 
in the expression for a^ neglects the synchronous rotation of 



Mj is the mass of Jupiter. 



the companion and its finite Love number. Including these 
effects leads to a slightly larger tidal shredding radius, im- 
plying that slightly more massive companions are required 
to unbind the CE and avoid tidal disruption. 

Since the stellar mass function is heavily weighted to- 
wards lower masses, we limit ourselves to a 1 progenitor. 
More massive primaries extend to larger radii (see Fig. [3| 
and swallow companions at farther distances, leading to 
wider period gaps. We consider two companions: a 1 Mj 
planet and a 10 A/j brown dwarf. 

For each evolutionary model and binary configuration, 
we calculate the minimum and maximum bounds of the sep- 
aration gap. Our results for rj — 0.7, 1 and 5, a; = and 1 
and Qo ~ 10^ and / = 1 are summarized in Table [2] The 
maximum of amin and minimum of Omax yield the minimum 
gap expected for a M* = IMq , Mc = IMj system. Note that 
for pedagogical purposes we include the results for mass-loss 
only (no tides) in Table [2] However, when determining min- 
imum period gaps, we only consider systems in which tides 
are acting. It should be stressed that even in the absence of 
tides, a minimum period gap exists. The functional depen- 
dence of tides on separation (Eq. 2) largely acts to change 
the location of the outer boundary of the gap. Depending 
on tidal prescription, this shifts the outer boundary of the 
by a factor of ~2-3 at most (Table [2|. 

From Table [2) we see that there should be a paucity 
of Jupiter-mass companions with periods <270 days around 
white dwarfs. Additionally from Table [2j we see that there 
should be a paucity of 10 Mj companions with periods be- 
tween 0.1 days (0.003 AU) and 380 days (0.75 AU). This is 
consistent with the tentative detection of a ~2 Mj planet 
in a >4 year (>2.75 AU) orbit around the white dwarf CD 
66 ( [Mullally et al.|2008[ [20091 ). Future surveys for low-mass 
companions around white dwarfs will have the ability to 
confirm or refute our prediction of a gap. Several efforts are 



either recently completed or are currently underway (Far 



ihi et al.||2006| [Tremblay and Ber geron|[2007[ [Hoard et al 



2007 



Farihi et al.|2008 Farihi|20'09 1. Once the samples are 



sufficiently large, observational identification of period gaps 
could help to constrain aspects of mass-loss and tidal theo- 
ries. 

For each M^, = IMq model, we can calculate the mini- 
mum mass companion able to unbind the envelope and sur- 
vive a CE phase. This occurs when the binding energy of 
envelope is at a minimum. Note that, in multiple planet 
systems, several close companions may incur a CE phase. In 
conjunction, multiple lower-mass companions could poten- 
tially supply the same energy as a single larger mass object. 
Although we do not investigate such a scenario, multiple 
companions may successfully expel the envelope such that 
all or a subset survive. For the 77 = 0.7 model, the binding 
energy of the envelope is minimized when past the peak ra- 
dius on the AGB. An 8.4 Mj mass companion is sufficient 
to unbind the envelope and survive a CE phase. For the 
77 = 1 model, the envelope binding energy is minimized at 
the peak of the RGB. For this system, a 9.8 Mj companion 
is sufficient to survive a CE phase. The — 5 model does not 
achieve a sufficient core mass to undergo a helium flash. As 
such, its maximal radial extent is less than the rj = 0.7 and 
rj — 1 models while the envelope binding energy is reduced. 
In this case, a ~7.8 Mj mass companion is sufficient to sur- 
vive a CE phase. For all of our models, most companions 
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Figure 6. Maximum initial semimajor axis, ai maxi that enters a 
common envelope for 1, 10 and 100 Mj companions with 1, 1.5, 
2, 2.5 and 3 Mq primary stars. The left column in each figure 
presents results for x = [Q'^ = Qo] while the right column 
presents results for x = I [Q'^ = Qo x (n/Ho)]. The evolutionary 
models of the primary star were computed with a Reimers mass- 
loss parameter of r; = 1. The left figure is in units of AU while 
the right figure is in units of the maximum radius on the AGB 
(see Table 1). 
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0.65 



Figure 8. The predicted period gaps for a 1 Mq progenitor with 
1 A/j (top) and 10 M,j (bottom) companions. The symbols repre- 
sent different Remiers rj values for the various tidal prescriptions 
listed in Table[2] For the 1 Mj system, no companion survives CE 
evolution. Thus, we predict a paucity of 1 Mj companions with 
periods <270 days. For the 10 Mq system, several companions 
survive CE evolution and are located in short-period orbits. The 
predicted period gap occurs between ~0.1 and 380 days. 



Figure 7. Maximum initial semimajor axis, a, that enters a 
common envelope for 1, 10 and 100 Mj companions with 1, 2 
and 3 Mq primary stars. The left column presents results for 
a; = [Q^ = Qo] while the right presents results for x = 1 
[Q* = Qo ^ (n/Ho)]. The color bar is in astronomical units. The 
evolutionary models of the primary star were computed with a 
Reimers mass-loss parameter oi rj = 0.7 (left figure) and jj = 5 
(right figure). In general, less mass-loss leads to tidal engulfment 
at a slightly larger semimajor axis. 



more massive than ~10 Afj will supply sufficient energy to 
unbind the envelope and survive the CEP at greater orbital 
separations]^ 



7 WD TRANSIT DETECTABILITY AND 
GALEX 

If a massive planetary companion or a brown dwarf is en- 
gulfed and survives the common envelope phase, it could 
be left in a close, short-period orbit around the resultant 
WD. In such a scenario, if the orbit were oriented edge-on, 
it is conceivable that an entirely new regime of exoplan- 
etary transiting observations would allow for detection of 
the companion. In contr ast to the observations pion eered by 

( |2000f | where 



Henry et al. ( 2000 1 and 



Charbonneau et al. 



transiting Jupiter-sized planets cause dips of ~1% in the 
light curves of main sequence solar type stars, a massive 
Jupiter transiting a WD could cause a total eclipse (a 100% 
dip!). How likely are such detections? In this context, rather 
than the oft-quoted a priori probability of /a that an ex- 
oplanet's orbit would be oriented so as to transit a main 



Note that even companions whose masses are several tens of 
Jupiter masses will still be within a few solar radii after surviving 
inspiral. 



^ See http:/ /exoplanet.eu' for an up-to-date census of the now 
more than 60 known transiting planets, and of all other known 
exoplanets. 
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" Companions which do not survive a common envelope phase are denoted by a — . 
^ The rj = 5 model does not achieve a sufficient core mass to undergo a helium flash. 
The minimum mass companion required to unbind the envelope and survive a CEP. 



Table 2. Period gaps for binary systems with different tidal formalisms (Q^ oc 11^; / = 1; no tides). The column quantities from left 
to right are Reimers mass-loss parameter rj, A/* the mass of the primary star, Ale the mass of the companion, Qo [Eq. x values 

such that Qi = Qo(n/no)^ [Eq.[TT], 

^^i,crit) the critical initial radius below/above which the companion is/isn't engulfed {see Fig. |T|| , 
Mcore,!, the core mass for a^nim Afcore,2, the core mass for Omax, Qmin the minimum separation (from CE evolution), Omax the maximum 
separation (from orbital evolution) , Pmin the minimum period (from CE evolution) , Pmax the maximum period (from orbital evolution) . 
Companions that do not survive common envelope evolution are denoted by — . Note that the time of engulfment for the above cases 
does not necessarily coincide with the minimum binding energy of the envelope during its evolution. Therefore, we include results for 
companions entering the giant star at the time of minimum envelope binding energy (-Ebind,min)- For graphical representations of this 
table see Fig. |8] 



sequence star, the probability of a total eclipse in a system 
where Rp > i?, is 



Porbit ' ~ 0.1 — yg— 



(18) 



where Rj is Jupiter's radius, indicating that roughly one 
out of every ~10 such post-CE systems might have an orbit 
that is oriented so that the companion eclipses the WD once 
per orbit. The duration of the eclipses will be of order 



W~:^~2xl0-^f^^ 
ZTva \ u/Rq 

times the length of the year (P), or 



(19) 



a/Rr. 



Ateciipsc = <5cciipscXP ~ (4 min) {Rp / Rj) ^ I ^^^^^ ^^^^ .(20) 

At any given time, the fraction of WDs being eclipsed by 
companions is of order 



/eclipse ^ fwD,p X Porbit X (5eclip 



(21) 



2 X lO'^fwD.p 



Rp/Rj 
cl/Rq 



where /wD,p is the fraction of WDs with short-period com- 
panions around them and Porbit, ^eclipse, Rp, and a are 
typical values of Porbit, ^eclipse, Rp, and a, respectively. 
The Galaxy Evolution Explorer (GALEX) is an ultraviolet 
space telescope capable of observing many hot WDs. Since 
GALEX time-stamps photons, it should be relatively easy 
to search the GALEX archive for transits. In order to detect 
an eclipsing system, a WD must be observed for a duration 
of at least Ateciipsc, or at two different epochs separated by 
at least AtecUpae such that it "winks out" and then reap- 
pears. Depending on the number of WDs that GALEX has 
observed (and on the duration of observations), a search of 
the GALEX archive might either detect some WD eclipsing 
binary systems - contingent on the distribution of radii and 
orbital semimajor axes of massive planets and brown dwarfs 
around WDs - or place constraints on the population of 
short-period companions to WDs. 
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T.l Progenitors of Cataclysmic Variables 

Via loss of orbital angular momentum, the short period sys- 
tems that survive a CEP may be brought into Roche contact 
with the WD. Such systems may have been detected as cata- 
clysmic variables (CV) but not recognized as post-planetary 
systems. If the donor is a massive planet with a rocky core, 
it may be possible for the core to survive the CV phase. This 
would result in a WD-|-rocky core short period binary. 



8 CONCLUSIONS 

By utilizing stellar evolution models from the ZAMS 
through the post-MS, we have followed the orbital dynamics 
of binary systems in which the companion is a 1 Mj planet, 
10 Afj brown dwarf, or 100 Mj low-mass MS star. Our evo- 
lutionary models incorporate a range of mass-loss rates and 
stellar masses. Dynamically, the orbital evolution is subject 
to mass-loss (which acts to increase the separation) and tidal 
torques (which act to decrease the separation). 

We employ two commonly used tidal prescriptions to 
investigate the differences during the post-MS. In particular, 
we find the following: 



• The tidal theory of Goldreich (|1963|) under the assump- 
tion that - 



10^ — 10^ leads to relatively weak tides, 
often capturing companions whose initial orbits are within 
~ 0.6 X -Rmax, where -Rmax is the maximum stellar radius 
for an evolutionary model. Companions with initial orbits 
larger than 0.6i?max but smaller than Rmax can escape cap- 
ture since their orbits expand due to mass loss of the pri- 
maries. 



• The tidal theory of Zahn (19661, under the assump- 
tion that f — 1 ( Verbunt and Phinney|1995 1, leads to com- 
paratively strong tidal torques, often capturing companions 
within ~(2 — 3) x JJmax- Nevertheless, even this is less than 
the values of ~5 — 7 x -Rmax that have been quoted in the 
literature jSoker|1996||Debes and Sigurdsson|2002[|Moe and| 



De Marcol2006 1 



• The ratio fli^max/J?*, where ai,max is the maximum ra- 
dius that is tidally engulfed, is not constant. 

• Future determination of the observational period gap 
for low-mass companions around WDs should help place 
constraints on the actual tidal theory acting in the post-MS. 

For each tidal theory, we determined the maximum sep- 
aration for which companions might be tidally engulfed (i.e. 
plunge into the primary star). These results serve as initial 
conditions for the onset of the common envelope phase for 
low-mass companions. Previous population synthesis predic- 
tions for post-AGB stars and PNe can be refined by incor- 
porating the methods outlined in this paper. 

For companions that incur a CE, under the assumption 
of maximum orbital energy deposition to the common enve- 
lope, we determined the maximum orbital radius at which 
a companion survives the interaction. By following the or- 
bital evolution of the closest companion that evades tidal 
engulfment, we predict a period gap for each binary system. 

For a binary system consisting of a 1 Mq primary with a 
1 Afj companion, we predict a paucity of Jupiter-mass com- 
panions with period below ~270 days. For a 1 Mq primary 
with a 10 Af,j companion, the gap occurs between ~0.1 and 



~380 days corresponding to ~0. 003-0. 75 AU. Note that our 
estimated gaps are conservative and are obtained by find- 
ing the minimum gap that might be expected for a range 
of mass-loss rates and a range of assumptions about tidal 
dissipation. It is unlikely that the true gaps would be nar- 
rower than the ranges quoted above, but they easily could be 
wider. As our knowledge of stellar evolution and tidal dissi- 
pation improves, so will our estimates of the ranges for these 
gaps. Finally, we note that the results of surveys searching 
for low mass companions to white dwarfs might help to con- 
strain theories of both stellar evolution and tides. 
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